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In quantum mechanics, useful experiments require multiple measurements performed on the iden- 
tically prepared physical objects composing experimental ensembles. Experimental systems also 
suffer from environmental interference, and one should not assume that all objects in the experi- 
mental ensemble suffer interference identically from a single, uncontrolled environment. Here we 
present a framework for treating multiple quantum environments and fluctuations affecting only 
subsets of the experimental ensemble. We also discuss a kinematic effect of indistinguishability not 
applicable to closed systems. As an application, we treat inefficient photon scattering as an open 
system. We also create a toy model for the environmental interference suffered by systems under- 
going Rabi oscillations, and we find that this kinematic effect may explain the puzzling Excitation 
Induced Dephasing generally measured in experiments. 
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I. INTRODUCTION 



It is because quantum theory makes probabilis- 
tic predictions that all quantitatively useful experi- 
ments require repeated measurements performed on 
identically prepared physical objects [l|. These iden- 
tically prepared physical objects constitute what we 
will call the experimental ensemble. And one can- 
not assume that all ensemble members, while treated 
identically by an experimenter, are also treated iden- 
tically by an uncontrolled environment. 

First we will explain what is meant by the exper- 
imental ensemble. Then to treat multiple quantum 
environments we will introduce a framework that 
one can use in addition to the typical reduced dy- 
namics of open systems 0. We will also explore a 
kinematic effect arising from the indistinguishabil- 
ity of different ensemble members. As applications, 
and to isolate the kinematics, we will work within 
the framework to model systems with either triv- 
ial or simple closed-system dynamics. We will treat 
as an open system the inefficient scattering of pho- 
tons from a beam splitter. Then we will create a 
toy model for real experimental systems undergoing 
Rabi oscillations. We find that the general yet of- 
ten puzzling phenomenon called Excitation Induced 
Dephasing results from the kinematics of indistin- 
guishable ensemble members, and consequently we 
find very good quantitative agreement with a wide 
variety of experiments. For one well-known experi- 
ment (H, we find unprecedented quantitative agree- 
ment. 



II. MEASUREMENTS AND THE 
EXPERIMENTAL ENSEMBLE 

The experimental ensemble is particularly easy 
to identify in modern experiments, such as 3-8], 
in which one repeatedly prepares, manipulates, and 
then measures the state of single physical objects, 
one at a time. 

For example, in Q, a single beryllium ion is con- 
fined in a Paul trap. It is prepared to be in a mo- 
tional and hyperfine ground state, and then a laser 
drives stimulated Raman transitions for a time du- 
ration, t, after which the state of the single ion is 
measured. Following that single observation, the 
preparation-evolution-measurement process is re- 
peated. For every value of t at which the experiment 
is performed, the result is the average of 4000 obser- 
vations [9(. There is an experimental ensemble for 
every measured value of t, and in this example every 
ensemble has size 4000. 

We offer these experiments only as clear illustra- 
tions of the experimental ensemble (though we will 
address Q and Q again in Section IVI) . Often mul- 
tiple objects from the experimental ensemble are 
present simultaneously in the laboratory. 

To generalize for systems and observables with 
discrete spectra, let A label a possible result of a 
quantum mechanical measurement. Assume that 
one performs N measurements, and that the result 
A is found n\ times. The results of these measure- 
ments are counting ratios: 



N ' 



(1) 



As is clear from experiments, the t in ([T]) corresponds 
to a duration in time rather than to a time on the 
laboratory's clock. Depending on the type of exper- 
iment, the ensemble size, N, may or may not be a 
function of t. 

The N measurements in (JT|) are performed on 
an experimental ensemble of N identically prepared 
physical objects. Because individual measurements 
are always assumed to be independent of each other, 
the individual members of the experimental ensem- 
ble are assumed not to interact with each other. 
Thus the experimental ensemble must not be con- 
fused with a many-particle state having interacting 
components, for which useful measurements require 
an experimental ensemble of many-particle objects. 
A major subject of this study is whether or not the 
members can be considered also statistically inde- 
pendent in terms of exchange symmetry or other 
correlations. This will be discussed in detail later. 



III. KINEMATIC ASPECTS OF QUANTUM 
THEORY 

One compares the experimental results in the form 
of ([1]) to theoretical probabilities. The theoretical 
image of any closed, quantum mechanical system 
is an operator algebra defined in a linear scalar- 
product space, $. The vectors <j> a span <£>, and every 
linear combination of the 4> a € 4> can possibly repre- 
sent the state of the system. Assume that prepara- 
ble states are represented by the basis vectors (p a . In 
practice, one infers from a measurement or from a 
preparation procedure the state of an experimental 
system. From the <pa one constructs the appropriate 
density operator, 

p(t) = ^w a \<p a (t))((j) a (t)\, with J~]w a = l, 



When comparing experimental results to theoret- 
ical predictions (|2J|, increasing N leads to stronger 
statistical conclusions regarding an inferred out- 
come, but increasing N does not change the outcome 
itself. In other words, the theoretical predictions on 
the right hand side of © are always independent of 
the number of members in an ensemble, N, because 



to represent the inferred state of the closed system. 
In this paper, we will only consider systems for which 
a takes discrete values. 



One characterizes an experimental observable by 
a prescription for how it is to be measured, and in 
the theory one represents it using a linear operator 
in $. For simplicity and because it is appropriate 
for the experiments described above, we will con- 
sider a projection operator, A = |A)(A|, into the one- 
dimensional subspace of $ that is associated with 
the observable labeled above by A. If A can be a 
value of the index a, then the subspace spanned by 
|A)(A| will be spanned also by one of the \(j> a ){4>a\, 
such that the inner product (X\(f) a ) oc 5 a \. It is 
straightforward to generalize to observables that do 
not have the <f> a as eigenvectors. 

These are kinematic aspects of quantum theory, 
and they are thus independent of any particular 
quantum system. The dynamics specific to a par- 
ticular system enters the theory only through a con- 
crete choice of the space $ and of the operator alge- 
bra defined in it. 



For projection operators such as A, the expecta- 
tion value is the probability of finding the discrete 
value A in a measurement. The predictions of quan- 
tum theory are calculated in the form of Born prob- 
abilities, Pa(p(£))j where 

V A {p(t))=Ti(Ap(t)). 

The explicit comparison between experiment and 
theory is 



(2) 



I 

no theoretical prediction should depend on the num- 
ber of times one performs an experiment to verify 
it. A theory's independence from N is a necessary 
condition for the identification of the experimental 
ensemble. 



a a | 

counting ratio = ^ (/?(£)) = Born probability. 

Experiment Theory 
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IV. MULTIPLE QUANTUM 
ENVIRONMENTS FOR A SINGLE 
EXPERIMENT 

Experimental systems cannot be perfectly iso- 
lated, and actual measurements always reveal the 
experimental signatures of environmental interfer- 
ence. In quantum theory there exist well-understood 
mechanisms to treat the interference between a 
quantum system and its environment [2] . The result 
is often a reduced dynamics, which can be used to 
calculate for an open system an increase in entropy 
or a loss of coherence. 

One has typically not considered the possibility 
that different ensemble members might suffer inter- 
ference from different environmental systems or from 
fluctuations within a single environmental system. 
That this possibility is relevant is especially obvious 
for the experiments described above 043 1 m which 
no two ensemble members are ever present simulta- 
neously in the laboratory. Of course, because quan- 
tum theory describes physical phenomena indepen- 
dently of how many ensemble members are present 
simultaneously, then it must be relevant for all ex- 
periments and any number of simultaneously present 
members. 

By ignoring the possibility that different ensem- 
ble members suffer interference differently, one im- 
plicitly assumes either that the uncontrolled envi- 
ronment treats them all identically or that experi- 
menters can prepare and control the state of the en- 
tire environment. In practice, this is either unlikely 
or impossible. 



A. Different environments from the beginning 

The simplest possibility for treating multiple 
quantum environments is to associate, for the en- 
tirety of the experiment, different ensemble mem- 
bers with different environmental systems. We will 
use m to label these different environmental systems, 
E m . The density operator representing the state 
of E m will be labeled PE m - The state of the sub- 
set of ensemble members suffering interference from 
the environmental system E m will be represented by 
p m . One must also label the time evolution operator, 
U m (t), defined so that as a function of the duration, 
t, after active preparation, the time evolution of the 
composite of system and environment is 

U m (t) 

l J m 



For a measurement performed at t, the states of the 
ensemble members are represented by the density 
operators 

Pm (t) = Tr Em (U m (t) Pm ® p Em Ul(t)) , (3) 

where Tr^ m denotes a partial trace over the degrees 
of freedom of the environmental system, E m . This is 
the typical reduced dynamics for open systems 0, 
but we have added the index m to label the dy- 
namics corresponding to the different environmental 
systems. 

The density operator, p(t), to be used in compar- 
ison with experimental data @ is then 

M 

P(t) = ^ a mPm(t)- (4) 
m=0 

Here M is the number of different environments 
deemed necessary for a model. So that the left hand 
side of (01 is a valid density operator, the real num- 
bers a m in Q must satisfy the convexity require- 
ments 

M 

a rn > and a m = I. 

m— 



B. Interference events 

With any given ensemble member remains 
for the entirety of the experiment in the subset cor- 
responding to a given environmental system. This 
will be problematic if there occur fluctuations in the 
environmental systems themselves, and the fluctua- 
tions affect only portions of subsets. We have as- 
sumed that members of the entire ensemble suffer 
interference differently, so we must assume also that 
members of any subset from (jj) suffer differently. In 
a general treatment, there will be no fixed subsets 
of ensemble members. 

To treat the possibility that members of any sub- 
set suffer environmental interference differently, one 
can assume that different fluctuations occur during 
different events throughout an experiment's dura- 
tion. Such events could describe something as sim- 
ple as a temporary loss of homogeneity of a magnetic 
field. Another possibility for an interference event is 
the actual passive measurement of an observable in 
the state of an ensemble member. In the context of 
decoherence, Leggett has called this measurement 
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process "garbling" 10] Q 

We imagine the following sequence of interference 
events. Recall that, when treating actual data ©, 
and thus when representing the state of ensemble 
members, the t represents a duration from active 
preparation. All members are therefore initially pre- 
pared at t = 0. 

1. A physicist actively prepares the initial state of 
all experimental ensemble members such that 
they are represented by a density operator la- 
beled po(t). At t = t* there occurs an inter- 
ference event that cannot be described by the 
closed-system dynamics. If we ignore for now 
any finite duration of the event and any other 
source of environmental interference, then be- 
fore the interference event (0 < t < £*) all 
ensemble members are represented by po(t). 

2. At t = t* the interference event affects some 
subset of the experimental ensemble. Because 
the closed-system dynamics cannot describe 
the interference, the affected members can no 
longer be represented by po(t). Thus, to rep- 
resent the state of the affected subset, one re- 
quires a new density operator: p\(t — t\). The 
unaffected members will continue to be in a 
state represented by po(t). 

3. If there is another interference event occurring 
at a duration ti,, then another density opera- 
tor, p2(t — t%), will be required. 

The density operator, p(t), to be used in compar- 
ison with experimental data ([2]) will thus exhibit a 
branching behavior as time progresses: 

fpoOo) 0<t<tl 
p(t) = < a p (t ) + oi pi(h) t\<i< t* 2 

[a po{t ) +axpi(h) +a 2 /9 2 (t2) t* 2 < t. 

(5) 

There is a notational complication here because the 
time parameter for experimental results, and thus 
the t in ^ , represents a duration from preparation. 
We will write t on the left hand side to represent 



Dirac has described the act of measurement as causing 
"a jump in the state of the dynamical system" While 
other interpretations exist, in active measurements ensem- 
ble members are never observed to be in superpositions. 
If nature does not distinguish between active and passive 
measurements, it is precisely this "jump in the state" that 
requires a new density operator. 



duration from active preparation of ensemble mem- 
bers. The branches on the right hand side of ([5]) will 
be functions of the different £j = t — i| , representing 
duration from the times t* when the pi became nec- 
essary. If the distinction is unnecessary, we will use 
t. For closed systems, this notation is not needed 
because t j — > t = t. 

Every line in ([5]) could be written in the same form 
as (|4]) , but here M would be the number of interfer- 
ence events. If there is a sequence of such events, say 
at t*, then the number of density operators on the 
right hand side will depend on the time duration, t, 
and more and more branches will appear as the time 
duration grows. 

Note that nothing new has been introduced to 
standard quantum mechanics, and that we have not 
assumed any new microphysical phenomena^ In 
the absence of environmental interference, these ap- 
proaches obviously reduce to the theory for closed 
systems. In the limit that all members suffer exter- 
nal interference identically, one has the typical for- 
malism for open systems. All we have done differ- 
ently is apply the standard theory to open systems in 
such a way that one can address the possibility that 
different ensemble members may suffer environmen- 
tal interference differently. One should not read ([5]) 
to mean that the density operator undergoes a par- 
ticular, dynamical change. Rather, this expression 
means that subsets of ensemble members suffer dif- 
ferently and therefore require different density oper- 
ators to represent their states. 

C. The branching process 

To represent the possibility that between t and 
t + At and before active measurement, some sub- 
set of the experimental ensemble suffers interference 
and afterward must be described by a new density 
operator, pi, let us write from ([5]) 

Po(t) — -> a aPa (t + At) +oi pi(t + At). (6) 

Note that © is not mathematically correct. Let 
Pi be constructed from vectors in the space In $i 
is defined the operator algebra representing the sys- 
tem's dynamics following one particular type of en- 



2 If the interference event is a passive measurement, the state 
vector collapse or wave function collapse [ill mechanism 
clearly cannot describe the branching process because it 
changes the density operator in {5j in ways that depend on 
the ensemble size, N. 
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vironmental interference, which may not be defined 
in the space of the closed system itself. The dy- 
namics of environmental interference is uncontrolled 
and may be unknown. In general, $^ ^ and 
one cannot perform the sum in ([6]). Instead, ([6]) is 
a phenomenological statement that is required be- 
cause one can neither perfectly isolate experimental 
systems nor ensure that all ensemble members suffer 
interference identically. 

While (JSJ) might not make mathematical sense, it 
will have a sensible physical interpretation if the cal- 
culated probabilities are real numbers between and 
1. The corresponding expression for the Born prob- 
ability is 

Tr(AopoW) — > a Tr(A p (i + At)) 

+ aiTr(A lPl (t + At)). (7) 

To calculate Born probabilities one must decide how 
the operators Aj or pi might be extended or limited 
to other spaces such that this requirement is fulfilled. 
For the remainder of this paper, we will assume that 
this can be done0 Again, the weights, do and ai, 
in ^ and ([7]) must satisfy ao + a\ = 1. 

We have ignored any dynamics hidden by the ar- 
row in ©. If the timescale associated with the in- 
terference event, At in ([5]), is very small compared 
to the dynamical timescales of the experiment, then 
one can safely ignore it. If it is not, then in models 
one must include interference dynamics along with 
the regular system dynamics. 



D. Partitioning and indistinguishability 

Finally, let us consider the weights, a*, in (@| 
and (|5|). The ai express the relative sizes of the 
ensemble's different subsets, which are in states rep- 
resented by the different pi. To determine the rel- 
ative sizes of these subsets, one must partition an 
ensemble of arbitrary size, N. 

To complicate matters, there is a kinematic re- 
quirement that limits one's knowledge of which 
members suffer interference from which environ- 
ments. When the spatial wave functions of n simul- 
taneously present ensemble members overlap, any 
n-body wave function must transform as one of the 



3 One possibility of interest is that Aj is the zero operator 
in In the beam splitter application we will demonstrate 
that this corresponds to inefficiencies in measurement or to 
dissipative effects. 



irreducible representations of the permutation group 
of order (nl). Even if the non-interacting ensem- 
ble members are sufficiently spaced that particle 
statistics can be ignored, when an external system 
cannot distinguish them, interactions still must re- 
spect a permutation symmetry. It is this condition 
that leads to the well-known phenomenon of super- 
radiance [12j. This phenomenon has previously been 
investigated in the context of open systems [13| and 
called collective decoherence, in which a bath (envi- 
ronment) cannot distinguish between n qubits. Here 
we have generalized to allow multiple baths, and we 
will attempt to treat the case that n varies and is 
not known. 

Clearly, if one does not know n, then one cannot 
explicitly address a permutation symmetry when de- 
termining the a.i . The physical interpretation of the 
permutation symmetry is that some simultaneously 
present members of the ensemble are indistinguish- 
able and cannot be given physically meaningful la- 
bels. 

We conclude with a set of restrictions on parti- 
tioning: 

1. p represents the state of an arbitrary number, 
N, of ensemble members (see Section IPTI]) . 

2. Subsets are therefore also of arbitrary size. 

3. One cannot explicitly address the number of 
ensemble members simultaneously affected by 
a single, uncontrolled environment. 

4. One cannot give individual ensemble members 
a physically meaningful label. 

These restrictions seem to prevent one from avoid- 
ing partitioning by defining within a standard re- 
duced dynamics any single, "effective" environment 
or any interaction operator that dynamically mimics 
the kinematic effects of indistinguishability for the 
average ensemble member. In the next two sections, 
we will use the framework to treat different physi- 
cal systems. To isolate the effects of kinematics, in 
Section [V] we create a model with trivial dynamics, 
and in Section I VII we create a model with simple 
dynamics. 

V. APPLICATION: BEAM SPLITTER 

Consider the simple system of a photon, 7, inci- 
dent on a beam splitter, being either transmitted 
into one phototube or else reflected into a different 
phototube. Figure [1] is a schematic of the experi- 
ment. The photon scattering system is interesting 
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be a\ and Sj, where i = 1,2, and 



® 



FIG. 1. Schematic of the photon scattering experiment. 
It is divided into System A, the incident photons, and 
System B, the reflected and transmitted photons. Re- 
flected photons are registered by phototube 1, and trans- 
mitted photons are registered by phototube 2. Not 
shown is an imagined source incident from the top. 



as the basis of a photon current. The probabilistic 
nature of the scattering at the bea m sp litter is for the 
photon current a source of noise [14| characterized 
by correlations between photons in the transmitted 
and reflected branches. 



A. Second quantization calculation 

Second quantization provides a simple and effec- 
tive framework for calculating correlations, in terms 
of occupation numbers in the incident, transmitted, 
and reflected currents. In the next section we will 
use our framework and compare to the results of 
this section. We are not interested in the micro- 
scopic scattering dynamics for a photon incident on 
a beam splitter, so following [l5[ for this simple case 
we will deal only with asymptotic in and out states 
at fixed energy, and ignore other quantum numbers. 

One comment is required, however. Because the 
scattering operator, S, does not change the space of 
asymptotic states, one requires another input chan- 
nel in Figure [IJ corresponding to photons incident 
from the top. Such photons, were they present, 
would transmit into phototube 1 and reflect into 
phototube 2. This imagined photon source is not 
present in the actual experiment, however, so one 
fixes its input as the vacuum state. Though it emits 
no photons, only when the imagined photon source 
is present in the calculation does one calculate the 
correct values for the correlation functions [15j . 

To preserve the symmetry implied by the theoret- 
ical approach, let us temporarily change the input 
system to include this imagined source. The creation 
and annihilation operators for incoming photons will 



\ai,a' A 



(8) 



The actual source, seen in Figure [TJ will correspond 
to i = 1 , and the imagined source will correspond to 
i = 2. 

For photons in the reflected and transmitted cur- 
rents, one defines the creation and annihilation op- 
erators, b\ and 6,, with i = 1,2 corresponding to 
the label on the phototubes in Figure [TJ and also 
satisfying flSJ). 

The reflection or transmission at the beam splitter 
is then treated as a scattering event, and one can 
relate the incoming photons to the transmitted and 
reflected photons by [14| . 



til r 2 2 



a 1 

a-2 



(9) 



Here the unitary scattering matrix is written in 
terms of the complex reflection coefficient, r, and 
the complex transmission coefficient, t, which sat- 
isfy |r| 2 + \t\ 2 = 1. The vacuum state for the input 
photons is |0)i|0)2 = |0). We will consider an in- 
put state, \n), containing n photons emitted by the 
actual source, i = 1, so that 



\n) 



|0). 



(10) 



Following [TJ], we let |r n | 2 = R and \t 21 \ 2 = T, and 



we assume unit efficiencies for measurement. The 
operator representing the occupation number in the 
reflected channel of Figure [JJ is hn = b\b\, and the 
the operator representing the occupation number in 
the transmitted channel of Figure Q] is ut = b\b2- 
Using © and (|10l) . one calculates for the correla- 
tions between occupation numbers: 

(hnfiT) — (n\b\biblb2\n) — RTn(n — 1) 
(fin) = nR 

(h T ) = nT (11) 
((Ah R ) 2 ) = ((An T ) 2 ) = nRT 
(An R An T ) = -nRT 



B. Scattering as an interference event 

In this section we will begin to work within the 
open systems framework of Section IIV1 and in the 
next section we will introduce losses. The scattering 
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system of Figure [T] is naturally divided at the beam 
splitter, into subsystems: 

1. System A, which includes the incident 
photons: 

The density operator pa (t) represents the state 
of the incident photons. It is defined in the lin- 
ear scalar-product space $,4, which is spanned 
by the ket In this basis, pa = The 
Hamiltonian operator is Ha, and System B is 
external, so we assume trivial dynamics with 
[p A ,H A ]=0. 

2. System B, which includes the reflected 
and transmitted photons: 

The density operator for a single photon, 
PB(t), is defined in the space $b, which is 
spanned by the orthogonal kets \r), represent- 
ing the state of reflected photons, and \t), 
representing the state of transmitted photons. 
The Hamiltonian operator is Hb, and the dy- 
namics is again trivial, with [ps, Hb] = 0. 

To measure the expectation values in , one ac- 
tively prepares n photons simultaneously in the inci- 
dent channel. The n photons then scatter from the 
beam splitter, and one measures the state of the n 
photons by counting the numbers registered in pho- 
totubes 1 and 2. If one performs N measurements, 
then one has identically prepared the state of n pho- 
tons a total of N times. 

Every incident photon is actively prepared in a 
state represented by p A = Eventually ev- 

ery photon reaches the beam splitter and is scat- 
tered, which is the interference event during which 
the uncontrolled environment passively prepares ev- 
ery photon such that its state must be described by 
ps- According to the framework proposed in Sec- 
tion llVl in addition to the typical reduced dynamics, 
density operators undergo the branching 

PA (t) w A p A (t + At)+w B p B (t + At) (12) 

whenever subsets of ensemble members are passively 
prepared by an uncontrolled environmento All dy- 
namics for in and out states is trivial, so we will drop 
the time dependence from equations in this section. 

The scattering matrix in ^ preserves photon 
number [lj|. From this we deduce for (fT2)) that 



wa + wb = 1. For the simple scattering experi- 
ment in which one assumes unit efficiencies, all pho- 
tons are eventually scattered into states in System 
B. We can therefore deduce further that wa = 0, 
fixing wb = 1- 

If one knows the scattering theory ([9]) for the pas- 
sive preparation of photons by the beam splitter, 
then one can write for n = 1 



PA^>PB = R\r)(r\+T\t)(t\ 

+r 11 t* 21 \r)(t\+t 21 r* 11 \t)(r\ 



(13) 



By assuming the most general form for ps, one 
would also arrive at ([T^]) up to the coefficients. 
For this example with trivial dynamics, the frame- 
work simply reproduces an expression (|13p com- 
monly used in the study of quantum optics [16j . 
The only difference for our framework is that (|I3[) 
represents the state of only the n = 1 subset of the 
experimental ensemble. 

Of course, one often will not know the complete 
theory for passive preparations by an unknown en- 
vironment. The measurements that lead to (jlip arc 
not sensitive to entanglement, and the dependence of 
the correlation functions on n is entirely kinematic. 

Without more sophisticated measurements or a 
complete scattering theory, by imagining that a 
physicist actively prepares any number of photons 
in a sequence nj, j = 1,2,3..., one can deduce a 
form for ps that, though incomplete, is consistent 
with ([IT]). For any sequence rij , the probability that 
any given photon is reflected is 



E 



R = Vr(pb) 



After an incident photon is scattered, then, the prob- 
ability to find it in the reflected channel is 



V R (pB)=Tr(\r)(r\p B ) = R. 



(14) 



Similarly, the probability to find the photon in the 
transmitted channel is 



V T (pB) = Tr(\t)(t\ P B)=T. 



(15) 



Knowing (j 14(1 and (|15p from measurements, and 
without knowing about the entanglement in ([To 
one deduces for (fl2ll 



PA —>PB = R\r)(r\+T\t)(t\. 



(16) 



4 Though pa and ps are not defined in the same space, for 
this scattering system with unit efficiencies, wa = 0, and 
we can ignore the complication for now. 



We are interested in exploring the kinematics of 
indistinguishability without knowing the number, n, 
of photons present, and therefore not concerned with 
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entanglement. If the beam splitter is to represent an 
unknown environment, when developing models one 
presumably would not even know about entangle- 
ment or the imagined photon source. For the re- 
mainder of this section we will ignore off-diagonal 
elements of the density operator (TIB")) and restrict 
ourselves to operators of the form (TT6")) . 

In (TIB|) we have also ignored the kinematics as- 
sociated with the simultaneous presence of n in- 
cident photons. One cannot start with (|16|) and 
directly calculate the correlation functions in (TlTj) . 
One way to include the kinematic information is to 
write the analogue of (|16[) for n incident photons. 
This is demonstrated in Appendix [X] for n = 3. 
Here we will show for general n that by counting 
indistinguishable reflection events and indistinguish- 
able transmission events, one can instead partition 
the experimental ensemble represented by p R in (1161) 
into an ensemble of n-photon objects while respect- 
ing permutation symmetry and the requirements of 
Section HVDl 

Each of the single photons is passively prepared 
with probability R to be in the state \r) and with 
probability T to be in the state \t). The ensem- 
ble members are otherwise indistinguishable, so to 



as expected from (|11|) . Though we have not yet in- 
troduced losses or dynamics, we have treated the 
scattering event as an interference event, and we 
have partitioned an experimental ensemble. 



partition the ensemble we will use the binomial dis- 
tribution, 

to average over the partitions containing combina- 
tions of k indistinguishable photons reflected with 
probability R and [n — k) indistinguishable photons 
transmitted with probability T. 

The expectation values (fin) and (tit) are then 
related to the first moment of the binomial distribu- 
tion. We calculate 

(nR)=Y J (t)R k T n - k k = nR, (17) 

k=0 ^ ' 

and 

(fir) = (f) R k T n - k (n -k) = nT. (18) 

fc=0 ^ ' 

In this method, we put k inside the summation 
to represent the number of photons reflected. We 
put (n — fc) inside the summation to represent the 
number of photons transmitted. Therefore (ti^ttt) 
is 



(19) 



(20) 
(21) 
(22) 

I 

C. Open systems and inefficiencies 

With the generalization of the theory of open sys- 
tems and the partitioning method of Section IV Bl it 
is straightforward to introduce losses into the photon 
scattering system. We will consider two sources of 
inefficiency. The first occurs if photons are scattered 
toward the phototubes of Figure [TJ but the photo- 



n 

' n 



(n R n T ) = J2 ( h )R k T n - k k{n - k) = RTn(n - 1). 

fc=0 



For the remaining relations, we find 

n 

((An R ) 2 } = ( 'I ) R k T n ~ k (k - nRf = nRT 



k=0 



((An T ) 2 ) =^2^jR k T n - k ((n-k)-nT) = nRT 
(An R Ah T ) = 'jh ( 1 f\R k T n - k (k - nR)((n - fc) - nT) = -nRT, 

k=Q ^ ' 
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tubes fail to register the photons. Unless one models 
the detector quantum mechanically, detector ineffi- 
ciency does not require an open system treatment. 
Let A' R be the operator representing observation of 
reflected photons, and let represent observation 
of transmitted photons. If there is some probability 
of a failed registration, then one has 

A' fl =>(l- efl )|r)(r| and (1 - e r )|t><t|, 

where cr is the probability that a photon fails to 
register in phototube 1, and ex is the probability 
that a photon fails to register in phototube 2. The 
probability to measure a reflected photon is 

V R ( PB ) = (l-e R )R, 

and the probability to measure a transmitted photon 
is 

V t {pb) = (l-e T )T. 

The second type of inefficiency occurs when the 
beam splitter scatters photons into such directions 
that they will not hit either phototube. Then one 
introduces another system, System C, that includes 
these photons. System B is the system of inter- 
est, and System C corresponds to the environmental 
system of the typical framework . Unmeasurable 
photons are in a state represented by a density oper- 
ator, pc, constructed from vectors in the space $c> 
which does not contain |r) or \t). Instead of p^|) . 
one uses 

p A ► W B PB + W C PC 

= w B (R\r)(r\+T\t)(t\) +w c pc (23) 

to represent passive preparation at the beam splitter 
(again ignoring entanglement). 

It is important to note again that (|23l) is not math- 
ematically correct because p B and pc are not defined 
in the same space. The probability after scattering 
to measure the observable represented by A is 

Vk b (wbPb) + V\ c . (w c pc) 
= w b V Ab (pb) + wcVk c (pc) , (24) 

where wb + wc — 1 , and Ap is the appropriate op- 
erator in the space $c- Because System C contains 
photons that are not measured, Ac is the zero oper- 
ator in $c". The probability to measure a reflected 



photon is 

Vr(wbPb) = w B R, 
and the probability to measure a transmitted photon 
is 

Vt{wbPb) = w B T. 

In real experiments, both sources of measurement 
inefficiency will appear, and one will have for the 
probabilities 

'Pr(pb) = w B (l - e R )R 

and 

Vt(pb) =wb(1- e T )T. 

Now we can check the results in (fTTj) , for the case 
that one does not have unit efficiencies. Let R' be 
the probability that a reflected photon is registered 
and T' be the probability that a transmitted pho- 
ton is registered. Let r\ be the probability that an 
ensemble member is not measured, such that 

R' + T' + i-/ = l. (25) 

If the experimental ensemble contains all photons 
actively prepared in the incident channel, then 

R' = w B (l-e R )R 

T' = w B (l-e T )T (26) 
rj = 1 - w B (l - e R R - e T T). 

If, instead, the experimental ensemble contains 
only those photons scattered into System B, then 
w B = I, and 

R' = (1 - e R )R 

T' = (1 - e T )T (27) 
r\ = €rR + e T T. 

In general, one will not know if a photon remains 
unmeasured because of inefficiencies in scattering or 
because of inefficiencies in the phototubes. In other 
words, one can measure R' and T", but one will not 
know which of ([2"B)) or (f2"T)) is appropriate. 

In (1231) we have the theory for the state of any sin- 
gle photon before and after an interference event. To 
calculate correlation functions for different n, again 
we partition the experimental ensemble. Now we 
must use the multinomial distribution: 
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n n n | 

= E E E ^w R,jT ' k ^ j = + r + 

j=0 k=0 £=0 ■>' ' ' 



n-X 



n n n | 

<^> - E E E ^w R,jT ' k ^ s ^+ k +* k = nT '( R ' + T '+v) 

9=0 fc=0 f=o ' ' 



71-1 



\n— 2 



j=0 fe=0 f=0 
n n n ■ 

(n«n T ) = E E E ^W R ' jT ' k ^ 6 -d+k+ijk = R'T'n(n - 1)(R> + T' + rj) 

j=0 k=0 1=0 

where 5 n .j+k+i is the Kronecker delta. Because of ()25|) . we have (fin) = nR' , (tit) — nT' , and {fiRjiT 
R'T'n(n - 1). 

Using the other expectation values are 

n n n | 

j=0 k=Q £=Q J ' 
n n n . 

^' — n ?,. n n J 



(28) 
(29) 
(30) 



n n n . 

j=0 fc=0 ^=0 J 
n n n | 

(An R An T ) = ^ E E "W^'^ C? - - = —nR'T' 

j=o k=o e=o 



!') 2 = nR'(T' + rj) 
") 2 - nT'(ii' + 77) 
- nR')(k - nT'] 



(31) 
(32) 
(33) 



Comparing ([28 ]) - (f33]) . where R' + T' + i] = 1, one 
can see how to generalize to imperfect measurements 
the results (HI]) of Section IVAl where R + T = 1 
and inefficiencies were ignored. The main difference 
is that ((An R ) 2 ) ^ ((Ah T ) 2 ) when rj ^ 0. If one 
writes from (jflj) that ((An fl ) 2 ) = riJ?(l - J?) and 
((An T ) 2 ) = nT(l - T), however, then one can cor- 
rectly take R -> R' and T -> T . 

In the language of the theory of open systems [|[ , 
the system of interest includes measurable photons 
in states described by ps , which is constructed from 
vectors in $3. The environment contains unmea- 
surable photons in states described by pc, which is 
constructed from vectors in $c- All scattered pho- 
tons are represented by vectors in $ = §b © 
which is the space of states of the combined "system 
plus environment," when entanglement with the en- 
vironment is ignored. 

In scattering experiments, one does not measure a 
dynamical dissipation in currents. In fp~6|) and ([23]) 
we therefore have models without dynamics. Ineffi- 
ciencies appear as a measurement event rate smaller 
than the incidence rate per area of incoming objects, 
for a given cross section. While an experimenter 
continues to supply photons, the losses appear only 
as measurement inefficiencies, which are treated by 
normalizing measured data to the total number of 
photons registered in the phototubes over the course 
of the experiment. 



VI. APPLICATION: RABI OSCILLATIONS 
EXPERIMENTS 

As an application with simple dynamics, we will 
create a toy model for systems undergoing Rabi os- 
cillations. Here time cannot be ignored, and we 
must deduce a form for environmental interactions. 
We will also implement the technique used in the 
previous section to partition the ensemble of single 
photons represented by (IT!)]) into an ensemble of n- 
photon objects. 



A. Rabi Oscillations 

The theory of Rabi oscillations is well estab- 
lished (l7| . It describes a two level system, with 
levels described by the projection operators \g)(g\ 
representing the ground state and |e)(e| representing 
the excited state. When the system is prepared at 
t = to be in the excited state, |e), and it is coupled 
to a correctly tuned radiation field, the probability 
at t to find the system in the ground state, \g), is 

TWlOK*)) - sin 2 (m) = l(l-oos(2fit)), (34) 

where fi is called the Rabi frequency. (The proba- 
bility to find the system in the same state in which 
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it was prepared is cos 2 (Qt).) These expressions are 
calculated using the dynamics for closed systems un- 
dergoing Rabi oscillations. The probability in (|3"1)) 
is to be compared with experimental results @- 

In actual experiments, because of environmental 
interference one never measures the result in . In 
this paper we wish to consider several very different 
but very clean and well-controlled experiments 0, [H, 

mm. 

In one experiment [3j, also discussed in Section [TT1 
internal levels of a 9 Be + ion couple to the harmonic 
binding potential. Rabi oscillations occur between 
two of the coupled internal and vibrational levels, 
and the oscillations between different sets of levels 
are measured. In another experiment Q, Rabi os- 
cillations are observed between the circular states of 
a Rydberg atom coupled to a field stored in a high 
Q cavity. In [l8j], the Rabi oscillations are between 
the spin states of two electrons in a double quantum 
dot. In [l9j], Rabi oscillations occur between mo- 
tional states of electrons bound to shallow donors 
in semiconductors. These motional states mimic the 
levels of a hydrogen atom. In three experiments (20I — 
HH, oscillations occur between excitation levels of 
electrons confined in quantum dots. 

The measured results in H H, EH , do not match 
the prediction in (I34[) . Instead, the measured prob- 
ability is fit by an appropriately damped sinusoid 
that is a function of the time duration, t, from ac- 
tive preparation, 

fWiM*)) = - e- 7t ~cos(2Q?)), (35) 

where 7 is an experimentally determined damping 
factor. 

In [l8j], amplitude, offset, and phase are fit as well, 
because the oscillations are measured also as a func- 
tion of a swept detuning voltage. In [22[ , pho- 
tocurrent is measured, and from it the dephasing of 
the Rabi oscillations is inferred. Numerically solv- 
ing the models they fit to their data, the subsystems 
undergoing Rabi oscillations fit (|35l . 



B. Dynamics for the toy model 

Already having the kinematic requirements for 
open systems and indistinguishable ensemble mem- 
bers, our next task is to deduce the proper dynam- 
ics for the experimental system in question. This 
means we must identify the spaces $j and the corre- 
sponding algebras of observables. Starting with the 
branched density operator from ([5]) , we calculate the 



probability at a duration t from active preparation: 

V A (p(t)) =Tr(A a Po {to)) + Tr(Ai ai pi(ii)) + . . . 

= a Q VA a (po(t )) +a 1 V Al ( Pl {t 1 )) + ..(36) 

where, as before, pi is constructed from vectors in 

For the long-time behavior of the measured result 
in (1331) . one has 

- lim V \9)(g\{p$)) = \- 

From this we deduce that there is no dissipation 
for the experimental system, and that the terms 
r PK i {pi{ti)) in (J2U) ought not to introduce the type 
of losses considered for the beam splitter in Sec- 
tion IV CI For these experiments, we then hypoth- 
esize that immediately after the interference event, 
$0 = $1 = for all i, so we can take <F t -> $. 
The algebra of observables will therefore be the typ- 
ical algebra for the closed system, so we will take 
Ai — > A = \g)(g\- We want a probability of the form 

M branch 

Pa(p(*))= a mV A (p m (t m )) (37) 

m=0 

where Mbranch is the number of branches. To keep 
track of the branches, we will find it convenient to 
retain the label i for the pt, even though all pi are 
defined in the same space, $. This formula is de- 
ceptively simple in appearance, because, in general, 
the number of branches, Mbranch, will depend on t, 
and the value of a m > for a fixed m! will depend on 

Mbranch • 

The pi are all defined in $, and to keep our model 
simple we will assume that during interference events 
the environment prepares the system to be in en- 
ergy eigenstates, \g) and |e), ignoring entanglement. 
When a passive preparation occurs at U = 0, the 
most general form for pi in ([5]) is then 

P i(U = 0) = Wi \g)(g\+Vi \e)(e\, (38) 

with Wi + Vi = 1. In (|3"5)) , w is a weight appropriate 
for ensemble members prepared passively to be in 
the ground state, and v is a weight for ensemble 
members prepared in the excited state. 

Note that this is not a new idea. Leggett has 
suggested [To| that decoherence "is exactly the re- 
sult of a 'measurement' whose result is uninspected," 
and he has called this process "garbling." With this 
toy model, developed within our framework, we will 
treat the effects of garbling. We ignore the dynamics 
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of the measurement process itself because it is sim- 
ply not understood and apparently occurs over very 
small durations in time. 



C. Born probability 

For simplicity, we will assume that the passive 
measurements for our systems are fair measure- 
ments, though this would not generally be required. 
To isolate kinematic effects, we will also ignore any 
reduced dynamics expressed in © and (@|. 

When modeling experiments or using more sophis- 
ticated models for environmental interaction, one 
would probably use a full simulation. To maintain 
transparency, however, we will search for an ana- 
lytical formula for our toy model. We will start 
with a sequence of interference events described in 
Section IIV Bl Because the damping of the oscilla- 
tions in is exponential, the decay rate is con- 
stant, and we can deduce that there should be a sin- 
gle timescale, At, for regularly spaced interference 
events. Then the Born probability when there may 
have been n interference events is V\ g )i g \ (p(nAt)) . 

As discussed in Section Mil models should not 
depend on the ensemble size or on the number of 
members simultaneously present. Intuitively, one 
would like a model that does depend on the num- 
ber of interference events that have occurred. Note 

I 



The explanation of (|39|) is straightforward. It is in 
the form of the probability in (|36|) , and each line cor- 
responds to a branch of the density operator in ([5])- 
From partitioning the ensemble for the beam splitter 
( Section [V}, we know that the size of a subset is pro- 
portional to the number of events in which ensemble 
members are appropriately prepared, if we count the 
events themselves as indistinguishable. We must re- 
late the binomial distribution to the passage of time, 
so at every step k, with < k < n — 4, we count the 



from our treatment of photon scattering that count- 
ing indistinguishable ensemble members is equiva- 
lent to counting as indistinguishable the events in 
which they are prepared. 

We are consequently led to a model with two pa- 
rameters. At is the timescale for interference. f3 rep- 
resents the probability that a randomly chosen time 
interval will have preceded an interference event, 
with < f3 < 1 and j3 = 1 for a perfectly isolated sys- 
tem. This can be better understood after equation 
(14"2")) . We will again use the binomial distribution, 

b(n,k,(})= (fjp k (l-p) n - k , 
and the normalization 

x:("V(i-/r- fe =i. 

fc=0 ^ ' 

To illustrate the model, let us write (155)) for the 
case that the average ensemble member will have 
suffered only one interference event before an active 
measurement occurs. This will truncate our formula 
at a reasonable size, and we will generalize to mul- 
tiple events below. The probability at 4At to find 
in \g) members that have been initially prepared at 
t = to be in |e) is 



+ 
+ 
+ 
+ 

(39) 

I 

(normalized) number of combinations for arranging 
the n time intervals between possible event times, 
such that k of them came before the single interfer- 
ence event. This number is given by 6(n, fc,/3), and 
is called in (|36j) and J5]). 

Because possible interference events occur at in- 
crements of the time scale, At, the weight 6(n, k, (3) 
is attached to any passive preparation occurring 
at kAt. Consider the second line, where k = 3, 
and compare to the notation in (|3"8|) . The weight 



V\ g){g \(p{4:At)) = 6(4, 4, /3) (sin 2 (ft 4At) cos 2 (ft At) + cos 2 (ft 4 At) sin 2 (ft At) 
6(4, 3, (3) (sin 2 (ft 3 At) cos 2 (ft lAt) + cos 2 (ft 3 At) sin 2 (ft lAt) 
6(4, 2, (3) (sin 2 (ft 2 At) cos 2 (ft 2 At) + cos 2 (ft 2 At) sin 2 (ft 2 At) 
6(4, 1, (3) (sin 2 (ft lAt) cos 2 (ft 3 At) + cos 2 (ft lAt) sin 2 (ft 3 At) 
6(4, 0, 13) (sin 2 (ft At) cos 2 (ft 4 At) + cos 2 (ft At) sin 2 (ft 4At) 
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w = sin 2 (r2 3At). The weight v = cos 2 (^3Af). The 
duration from passive preparation to AAt is then 
(4 — 3) At = lAt, so the functions with argument 
lAt contain the time dependence of the Born prob- 
ability of the mixed state. 



Let us introduce the notation V^l, g ,{p(nAtj) to 
represent the Born probability under the assumption 
that members on average will have suffered i inter- 
ference events before active measurement. Then for 
general n, 



V \l)(g\ (p( nA fy = & ( n ' fc ' i 8 ) ( sin2 (° kA t) cos 2 (fl(n - k)At) + cos 2 {Q kAt) sin 2 (^(n - k) At)) . (40) 

k=0 

By simply exchanging cos 2 (f2fcAt) and sin 2 (fi kAt), we calculate P| V,, (p(nAt)) , which is the probability 
to find the ensemble members in the excited state. 

In (I40[) , we have assumed that ensemble members will have suffered one interference event. To allow for the 
possibility of multiple events, the terms sin 2 (fife Af) and cos 2 (fl k At) must be replaced with new functions 
of kAt, that predict the effects of interference events prior to the single event assumed in (I4U1) . For general 
i, we get the recursive expression 

n 

V§ {g] (p(nAt)) =J2Kn,k^)(v^(p(kAt)) cos 2 ^-^^,^ (p(kAt)) sin 2 (tl(n-k) At)) . (41) 

k=0 



We have written the Born probability as a func- 
tion of nAt. The final step is to scale our result 
back to the continuous t. The first moment of the 
binomial distribution is 



(42) 



events, is therefore 



k=0 



After stepping through time to nAt, on average /3 n 
of the intervals will have preceded the interference 
event number i. This provides us a physical inter- 
pretation of our two parameters. To ensure that the 
time scales with something physical, as indicated by 
the exponential nature of the damping, and scales 
as well with the recursive index i, we will need to 
use (k)Ai = (3 nAt — t. After a calculation of the 
probability as a function of nAt, we make the re- 
placement 



t 



P At 



(43) 



This restricts us to non-zero values of (3 and At. We 
have also simply interpolated between the discrete 
values of n at which (j4T|) is actually defined. Because 
our time scale is understood to be an average value, 
it would be inappropriate to assume for our model 
anything more complicated. 

The probability at the duration, t, from active 
preparation and assuming on average i interference 



V 



(0 

Is) (si' 



Is) (si 



P At 



(44) 



In Figure [2] we have plotted "P^gh g i (p(t))- F° r clean 
experiments we expect a good approximation for i = 
5. The agreement with the experimentally measured 
damped sinusoid is very good. 



D. Measurable consequence of 
indistinguishability 



Experiments reveal that the damping factor, 7, 
seems to depend generally on the Rabi frequency, f2, 
regardless of the nature of the experiment. In the 
experiment with 9 Be + ions the different levels 
are described by the kets | i, n) and | f, n + 1) , where 
I l) and I f) are internal states of the Be ion, and \n) 
represents vibrational Fock states. Rabi oscillations 
are measured for the frequencies 0, Q 



202 p -0.202 2 /2 

0„, n+1 = ET-^ e ■ i r \(0.202 2 ), 

Vn + 1 



(45) 



where L„ is the generalized Laguerre polynomial. 
The corresponding damping factor, j n , is measured 
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5 10 15 20 25 30 
Dimensionless time coordinate, D.t 

FIG. 2. Plot of the Born probability, r P\gl g Ap(t)), for 
indistinguishable ensemble members. For the dots we 
have used equations (|4ip and (|44[1 . The solid line is a plot 
of the damped sinusoid (|35[) that fits the experimental 
data. We have used QAt ~ 0.7 and f3 — 0.995, and we 
have fit j/n = 0.039. 



to increase with n according to 

^«(l + n) - 7 . (Measured) (46) 
7o 

In the experiment [l8[ with spin states of two elec- 
trons in a double quantum dot, the damping factor, 
7, is stated to be proportional to the Rabi frequency, 
f2, though no mathematical relation is given. 

For pJrHH, this dependence has been called Ex- 
citation Induced Dephasing (EID), and its cause has 
remained an open question. In [2lL l22j | , the depen- 
dence of the damping factor on the Rabi frequency 
has been found to be 7 oc f2 2 . 

These results for such different experiments sug- 
gest that some dependence of 7 on may be due to 
a general, kinematic effect rather than any dynam- 
ics specific to an experimental system or its uncon- 
trolled environment. 

When we include the kinematics of indistinguisha- 
bility in our toy model, we find that the damping 
factor, 7, does in general depend on the Rabi fre- 
quency, f2, in agreement with measurements and the 
phenomenon called Excitation Induced Dephasing. 
For (3 close to 1, the probability is dominated by the 
terms proportional to b(n,k = n,/3). Solving the 



rather crude, truncated forrr|f| 

a 

V lg){gl (p(nAt)) « ^6(n,fc,/3)sin 2 (r!fcAf) (47) 

fc=0 

and using (|4"4")) . we get 

*Wi (/>(*)) * 1( 2 - - P{1 ^ 2lA " 2 ))^ 

-(1 - 0(1 - e +2lAin )) flfe) .(48) 
For small At, (ggj) is 

V\ g )(9\ (Pity = \ (l - e" 7 ' (cos(2m) + O(AP))) , 

(49) 

where 

7 = 2(l-0)fi 2 At + O(Af 3 ). (50) 

For j3 w 1 and ilAt -C 1, 7 is quadratic in fi. This 
matches the results in [2l|, |22j . 

Because we made no assumption regarding the 
dynamics of environmental interference, the same 
model is general enough to apply also to the sys- 
tem in Q, and with it we fit the measured relation 
(|46|) between 7 and f2. Figure [3] is the result of fit- 
ting the damping factor, 7„, to "P\g\/ g i (/°(*)) j calcu- 
lated with the sequence of frequencies in (|4"5|) and 
with il At ps 0.2. (The exponent can be shifted by 
choosing different time scales.) Unfortunately, we 
are limited by computational resources to n < 6, 
and for n = 6 we were only able to use the first few 
oscillations when fitting. 

Figure [3] shows a very good quantitative agree- 
ment with experiment (|46[) . Our study indicates 
that the decoherence measured in Q is explained 
simply by the environmental interference character- 
ized in Section IVIBI and having a characteristic fre- 
quency of « 5Qq- 

In the standard decoherence program, one typi- 
cally treats these experiments using the decoherence 
master equation, and ignores the possibility that 
different ensemble members may suffer environmen- 
tal interference differently. The generic solution of 
the master equation describing a system undergo- 
ing Rabi oscillations @ predicts a constant 7, with 
no dependence on ft. This is one reason why the 
measured results have been puzzling. 



This would correspond also to a model for passive prepa- 
rations always into the ground state. 
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Level Index, n 



FIG. 3. Matching the experimental results for the ratio 
of damping factors, ^ . The large dots result from our 
theoretical calculation of the Born probability using (|41|l 
and (144[) . The solid line is the experimentally measured 
relation, (1 + n) 0,7 . To indicate a scale for the exponent, 
the thin dashed line is a plot of (1 + n)°' a , and the thin 
dotted line is a plot of (1 + n) 0,6 . 



When assuming that all ensemble members suffer 
interference identically, one has not needed to con- 
sider the kinematics associated with indistinguisha- 
bility. As a consequence, one has searched for dy- 
namical explanations for the measured dependence 
of 7 on n. For example, for the experiment on Be 
ions Q, one has required models with detailed in- 
teraction terms exclusive to the experiment itself. 
Among other things (23l-l25|. one has assumed laser 
amplitude noise [26| , decay from intermediate levels 
and, simultaneously, char ge- c oupling of the ion with 
the trap's electrodes [27ll28j, or feedback damping 
from polarized background gas polarized by the os- 
cillating ions themselves [29|. While plausible, none 
of these studies has resulted in an agreement as 
quantitatively good as that in Figure |3l 

Because they contain dynamics exclusive to exper- 
iments on ions in a Paul trap, many of these models 
have in effect been tuned to that experiment. They 
are thus not applicable to the other types of exper- 
iments [ij-lll]. The models [HI, [H, [U (U used 
to explain the EID measured in experiments with 
shallow donors and quantum dots have typically as- 
sumed a feedback dynamics from, for example, ex- 
citation induced phonons in substrates, and these 
models have often been quantitatively successful. 

When one invokes dynamics to explain a kine- 
matic effect, of course, one risks overly tuning mod- 
els and misdiagnosing the underlying physics. We 
have demonstrated that none of these dynamical 
mechanisms is necessary to explain Excitation In- 



duced Dephasing. As shown in Appendix [B] an 
analogous model for distinguishable ensemble mem- 
bers predicts no dependence of 7 on tt. We there- 
fore conclude that a dependence of 7 on SI, or EID, 
may indeed be general, as suggested by experiment, 
and that it may be a kinematic effect of the indis- 
tinguishability of separate, uncontrolled interactions 
between quantum systems and their environment. 



VII. CONCLUSION 

In addition to the typical reduced dynamics for 
open quantum systems [2|, to understand environ- 
mental interference one must faithfully represent the 
physics involved in the preparation, evolution, and 
measurement of all members of the experimental en- 
semble. This means one must address the possibility 
that different subsets of the experimental ensemble 
may suffer environmental interference differently. 

We have introduced a framework to treat multiple 
quantum environments. We have discussed a kine- 
matic effect of indistinguishability that affects only 
open quantum systems. The indistinguishability of 
experimental ensemble members greatly complicates 
the partitioning that our approach requires. Based 
on the beam splitter model, we have made the sim- 
plifying hypothesis that environmental interference 
occurs in events that are themselves indistinguish- 
able. We have not yet investigated alternative treat- 
ments, but the kinematic nature of the effect would 
seem to preclude the possibility of treating it dy- 
namically through a single, "effective" environmen- 
tal system and a typical reduced dynamics. 

To demonstrate our generalization for multiple 
quantum environments, within the framework we 
have described two simple systems. We calculated 
correlation functions for photons scattered at a beam 
splitter, with and without losses. We created a de- 
coherence model for systems undergoing Rabi os- 
cillations. In one case Q, we have found unprece- 
dented quantitative agreement with measurements. 
We have also discovered that the kinematic effect of 
indistinguishability can explain the generally mea- 
sured Excitation Induced Dephasing that has previ- 
ously required different dynamical explanations for 
different experiments. Full treatment of experiments 
will likely require more detailed models address- 
ing multiple sources of environmental interference, 
dressed states, reduced dynamics, etc. But the qual- 
itative and quantitative success of our single, general 
model for different types of Rabi oscillations exper- 
iments is very promising. 
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Appendix A: Including kinematics directly for 
photon scattering 



Here we treat photon scattering as an open sys- 
tem and directly include the kinematics of indistin- 
guishability for a fixed value of n. This corresponds 
to partitioning an ensemble of an arbitrary number 
of single photons into an ensemble of an arbitrary 
number of n-photon objects. 



Let pA, n represent the state of the n actively pre- 
pared, incident photons. The density operator pA.n 
is defined in the tensor product space <&® n , which is 
spanned by the ket | i) 2 ■ ■ ■ \i)n = 1*1*2 ■ ■ ■ in 

). The 

operator pA.n therefore remains trivial even when n 
photons are incident. 



The beam splitter passively prepares every mem- 
ber of the experimental ensemble to be in a state rep- 
resented by pB,n, which is defined in the tensor prod- 
uct space Again we will ignore entanglement, 
so pb,u is a sum of projection operators into each of 
the 2 n one-dimensional subspaces of $® n . Each pro- 
jection operator is weighted by R nR T nT , where ur 
and riT are the numbers of reflected and transmitted 
photons, respectively, represented by the associated 
one-dimensional subspace. 



Consider, for example, the case where n — 3. One 
writes the properly symmetrized density operator, 
pB, 3, to be used on the right hand side of (|T2"j). as 



Pb,3 = R 3 \rrr)(rrr\ 

rrt)(rrt\ + \rtr)(rtr\ + \trr)(tr 

+RT 2 (\rtt)(rtt\ + \trt)(trt\ + \ttr)(ttr\ 
+T 3 \ttt)(ttt\. (Al) 
If R + T = 1, then Trp B , n = 1 for all n. 



In the same basis, it is straightforward to con- 
struct operators corresponding to the correlation 



functions calculated in Section IV Al For example, 
nn = 3 \rrr)(rrr\ 

+2\\rrt){rrt\ + \rtr)(rtr\ + \trr)(trr\j 

+l(\rtt}{rtt\ + \trt)(trt\ + \ttr}{ttr\j .(A2) 

Proceeding in this fashion, one can reproduce 
for n = 3. 

Appendix B: The Rabi oscillations model for 
distinguishable objects 

Here we develop a model based on ([5]) but for 
distinguishable ensemble members and interference 
events. Our scenario is as follows: 

1 . The members of the ensemble can again suffer 
environmental interactions at the times nAt, 
where n = 1,2,3... 

2. At every nAt, there is some probability, (1 — 
77), for a member to suffer perturbation and 
thereby to be prepared passively. 

We have again used At for the timescale of inter- 
ference. But now we have used the parameter, 77, 
with < 77 < 1, to represent the susceptibility of en- 
semble members to environmental interference. The 
parameter 77 is the probability for any given mem- 
ber of the ensemble not to suffer interference at one 
of the times nAt. For a perfectly isolated system, 
77=1. 

We will again assume that passive measurements 
are fair measurements, and we will not assume any 
reduced dynamics. We can then write a general for- 
mula for the Born probability as more and more 
branches form: 



(*>(*)) = 



p (t) 0<t<lAt 

lAt<t<2At 



p n {i) nAt <i < (n+ I) At 

( B1 ) 

will have the form of the sum 



Here, every p n (t) 
in ([37]), and M branch = n. 

If members are actively prepared to be in the ex- 
cited state, and because none will have suffered en- 
vironmental interference before t = 1 At, we have 
for the initial value po(t) = sin 2 (fit). It is straight- 
forward to find from the theory of Rabi oscillations 
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Pi{t) = VPo(t) + (1 - r7)(po(lAt)cos 2 (f7 (t- lAt)) + (l -p (lAi)) sin 2 (fl (t- LAt))). (B2) 
For general u, we have 

Pn(t) = VPn-ify + (1 - 77) (p„_ 1 (nAi) cos 2 (ft (t - nAt)) + (1 - p„_!(nAf)) sin 2 (ft (t - nAt))\ (B3) 



Note the reappearance of a recursive structure sim- 
ilar to that of (|4T|) . Recursively partitioning the 
experimental ensemble in this way (|B3[) . however, 
requires that one can at any time label different en- 
semble members and know to which partition they 
belong. Because the recursion index, n, is also the 
index counting time in functions of nAt, there is 
no need to rescale the time parameter, as was done 
in (021). 

Figure [4] shows the results of two sample calcula- 
tions using (|B1[) with (|B3[) . The dots are calculated 
from our model. The solid lines are plots of the 



decaying sinusoid in (|35[) , which fits the experimen- 
tal measurements. For both figur es we have used 
flAi ss 0.08. The results in Figure 4(a) were calcu- 
lated using 77 = 0.99 and resulted in a fitted value 
for the damping factor of j/fl = 0.05. The results 
in Figure [4(b) were calculated using r\ = 0.997 and 



resulted in a fitted value for the damping factor of 
7/O = 0.015. Recall that r\ = 1 for a perfectly iso- 
lated system. 

This calculation for distinguishable members re- 
sults in a fitted damping factor, 7, that is indepen- 
dent of 17. 
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